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We investigate the mass spectrum of the pentaquark baryon (ududs) in quenched 
lattice QCD with exact chiral symmetry. Using 3 different interpolating operators, 
we measure their 3x3 correlation matrix and obtain the eigenvalues A+ (t) with ± 
parity, for 100 gauge configurations generated with Wilson gauge action at f) = 6.1 
on the 20 3 X 40 lattice. For the lowest-lying J p = l/2 — state, its effective mass is 
almost identical to that of the KN s-wave, while for the lowest-lying J p = 1/2+ 
state, its effective mass is smaller than that of the KN p-wave, especially for the 
regime m u < m s . By chiral extrapolation (linear in m^.) to = 135 MeV, 
we obatin the masses of the lowest-lying states: m(l/2 — ) = 1424(57) MeV and 
m(l/2 + ) = 1562(121) MeV, in agreement with the masses of mjf + mpf — 1430 
MeV and 6+(1540) respectively. 



1. Introduction 

The recent experimental observation of the exotic baryon + (154O) (with 
the quantum numbers of K + n) by LEPS collaboration 1 at Spring-8 and the 
subsequent confirmation from some experimental groups has become one of 
the most interesting topics in hadron physics. The remarkable features of 
9 + (1540) are its strangeness S = +1, and its exceptionally narrow decay 
width (< 15 MeV) even though it is ~ 100 MeV above the KN threshold. 
Its strangeness S = +1 immediately implies that it cannot be an ordinary 
baryon composed of three quarks. Its minimal quark content is ududs. 
Nevertheless, there are quite a number of experiments 2 which so far have 
not observed O + (1540) or any pentaquarks. This casts some doubts about 
the existence of 6+ (1540). 

Historically, the experimental search for O + (1540) was motivated by 
the predictions of the chiral-soliton model 3 , an outgrowth of the Skyrme 



"This work was supported in part by the National Science Council, ROC. 



1 



February 2, 2008 2:4 Proceedings Trim Size: 9in x 6in 



5Q04 



2 

model. Even though the chiral solition model seems to provide very close 
predictions for the mass and the width of 9+ (1540), obviously, it cannot 
reproduce all aspects of QCD, the fundamental theory of strong interac- 
tions. Now the central theoretical question is whether the spectrum of QCD 
possesses + with the correct quantum numbers, mass, and decay width. 

At present, the most viable approach to solve QCD nonperturbatively 
from the first principles is lattice QCD. Explicitly, one needs to construct an 
interpolating operator which has a significant overlap with the pentaquark 
baryon states. Then one computes the time-correlation function of this 
interpolating operator, and from which to extract the masses of its even and 
odd parity states respectively. However, since any (ududs) operator must 
couple to hadronic states with the same quantum numbers, it is necessary 
to disentangle the lowest lying pentaquark states from the KN scattering 
states, as well as the excited pentaquark states. 

To this end, we adopt the variational method 18 ' 19 , and use three differ- 
ent interpolating operators of (ududs) to compute their 3x3 correlation 
matrix, and from its eigenvalues to extract the masses of the even and 
odd parity states. These three interpolating operators (with 1 = 0) are 

4,5,6,7,8,9,10. 

-S a; /3e(75)/3 t; d a ; ne (75u) a ; QC } (1) 

(.Oz)xa = [u T Cj 5 d] xc {s Xl g e (7 5 ) / 3 7) u x?)C (7 5 d) :i;ae 

(75 u) 

} (2) 

(Os) xa = e cde [u T C l5 d] xc [u T Cd] xd (Cs T ) xae (3) 

where u, d and s denote the quark fields; e a tc is the completely antisym- 
metric tensor; x, {a, b, c} and {a, 13, 77} denote the lattice site, color, and 
Dirac indices respectively; and C is the charge conjugation operator. Here 
the diquark operator is defined as 

[u T rd]:m = e a bc(u xa bl 'apdxfSc ~ d xab T a fjU x fj c ) (4) 

where T a p = — Tp a . Thus the diquark transforms like a spin singlet (l s ), 
color anti-triplet (3 C ), and flavor anti-triplet (3/). For Y = C75, it trans- 
forms as a scalar, while for F = C, it transforms like a pscudoscalar. Note 
that 0\ , O2 , and O3 all transform like an even operator under parity. 

In the Jaffe-Wilzcek model 7 , each pair of [ud] form a diquark. Then 
the pentaquark baryon 9 ( [ud] [ud] s) emerges as the color singlet in (3 C x 
3 C ) x 3 C = l c + 8 C + 8 C + 10 c , and a member (with S = +1 and / = 0) 
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of the flavor anti-decuplet in 3/ x 3/ x 3/ = If + 8/ + 8/ + 10/. Now, 
if one attempts to construct a local interpolating operator for [ttci][ud]s, 
then these two identical diquarks must be chosen to transform differently 
(i.e., one scalar and one pseudoscalar), otherwise e a b c [ud\ x b[ud\ xc s xaa is 
identically zero since diquarks arc bosons. Thus, when the orbital angular 
momentum of this scalar-pseudoscalar-antifermion system is zero (i.e., the 
lowest lying state) , its parity is even rather than odd. Alternatively, if these 
two diquarks are located at two different sites, then both diquark operators 
can be chosen to be scalar, however, they must be antisymmetric in space, 
i.e., with odd integer orbital angular momentum. Thus the parity of lowest 
lying state of this scalar-scalar-antifermion system is even, as suggested in 
the original Jaffe-Wilzcek model. (Note that all correlated quark models 
e.g., Karliner-Lipkin model 11 and flavor-spin model 12 , advocate that the 
parity of 9 + (f540) is positive.) 



2. Computation of quark propagators 

Now it is straightforward to work out the pentaquark propagator (® X a®ys) 
in terms of quark propagators. In lattice QCD with exact chiral symmetry, 
quark propagator 13 with bare mass m q is of the form (D c + m 9 ) _1 where 
D c is exactly chirally symmetric at finite lattice spacing. In the continuum 
limit, (D c + rriq)^ 1 reproduces [j^(d^ + iA^) + rriq]^ 1 . For the optimal 
domain-wall fermion 14 with N s + 2 sites in the fifth dimension, 

l + l5 S(H w ) ofII i-n^ir. 



D c = 2m S(H W )- 

_ I - LO s H w _ 
1 s 1 „ , tl w — "f5-U w , 

1 + Lu s ti w 

where D w is the standard Wilson Dirac operator plus a negative parameter 
— m (0 < m < 2), and {uj s } are a set of weights specified by an exact 
formula such that D c possesses the optimal chiral symmetry 14 . Since 

(D c + m,)- 1 = (1 - rm^iD-^mg) - r], r = ^~ 

and D{m q ) = m q + (mo — m q /2)[l + -f5S(H w )], thus the quark propagator 
can be obtained by solving the system D(m q )Y — I with nested conjugate 
gradient 15 , which turns out to be highly efficient (in terms of the precision 
of chirality versus CPU time and memory storage) if the inner conjugate 
gradient loop is iterated with Neuberger's double pass algorithm 16 . 
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We generate 100 gauge configurations with Wilson gauge action at (3 = 
6.1 on the 20 3 x 40 lattice. Then we compute two sets of (point-to-point) 
quark propagators, for periodic and antiperiodic boundary conditions in 
the time direction respectively. Here the boundary condition in any spatial 
direction is always periodic. Now we use the averaged quark propagator 
to compute the time correlation function of any hadronic observable such 
that the effects due to finite T can be largely reduced 17 . 

Fixing mo = 1.3, we project out 16 low- lying eigenmodes of \H W \ and 
perform the nested conjugate gradient in the complement of the vector 
space spanned by these eigenmodes. For N s — 128, the weights {iv s } are 
fixed with A mm = 0.18 and \ max = 6.3, where A mm < X(\H W \) < \ max 
for all gauge configurations. For each configuration, (point to point) quark 
propagators are computed for 30 bare quark masses in the range 0.03 < 
m q a < 0.8, with stopping criteria 10~ n and 2 x 10~ 12 for the outer and 
inner conjugate gradient loops respectively. 

Then the norm of the residual vector of each column of the quark 
propagator is less than 2 x 10~ n , \\{D C + m q )Y — I|| < 2 x 10 -11 , and 
the chiral symmetry breaking due to finite N s (— 128) is less than 10~ 14 , 
a = \ S 2 Y/Y^Y — l| < 10~ 14 , for every iteration of the nested conjugate 
gradient. 

3. Determination of a -1 and m s 

After the quark propagators have been computed, we first measure the 
pion propagator and its time correlation function, and extract the pion 
mass (m^a) and the pion decay constant (f n a). With the experimental 
input U = 132 MeV, we determine a" 1 = 2.237(76) GeV. 

The bare mass of strange quark is determined by extracting the mass 
of vector meson from the time correlation function 



At m q a = 0.08, M v a = 0.4601(44), which gives My = 1029(10) MeV, in 
good agreement with the mass of 0(1020). Thus we take the strange quark 
bare mass to be m s a — 0.08. Then we have 10 quark masses smaller than 
m s , i.e., m u a = 0.03,0.035,0.04,0.045,0.05,0.055,0.06,0.065,0.07,0.075. 
In this paper, we work in the isospin limit m u = md- 




(U=l s 
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4. The 3x3 correlation matrix for 

Next we compute the propagators ((Oi) xa (Oj) y s) with fixed y — (0, 0), and 
their time correlation functions Cfj (t) with ± parity 



cfAt) 



E 



tr 



1±74 



(0,(^)0,(0,0))/ 



where the trace sums over the Dirac space, and the subscripts / and U 
denote fermionic average and gauge field ensemble average respectively. 
Then the 3 x 3 correlation matrix C ± (t) = {Cy(*)} can be constructed. 
Now with a judiciously chosen to, we diagonalize the normalized correlation 
matrix C ± (*o) _1/2 ± (i)C ± (i )~ 1/2 and obtain eigenvalues {Af(t)}, and 
from which to extract the masses {mf } of the lowest lying and two excited 
states for ± parity respectively. This is the variational method 18,19 to 
disentangle the lowest and the excited states. Then the mass mf can be 
extracted by single exponential fit to Af(t), for the range of t in which the 
effective mass M Q g(t) = \n[A(t)/A{t + 1)] attains a plateau. 



p 1 + 

Lowest lying state J = j 



m„a = 0.1 




Effective mass of + 
lowest lying state J p = ^ 



10 

t 



15 20 



(a) 



(b) 



Figure 1. (a) The eigenvalue A+(t) of the lowest positive parity state, for m u a = 0.1. 
The solid line is the single exponential fit for 9 < t < 14. (b) The effective mass 
M eS (t) = \n[A{t)/A(t + 1)] of A+(t) in Fig. la. 



In Fig. la, the eigenvalue A + (t) corresponding to the lowest J p = 1/2+ 
state is plotted versus the time slices, for m u a = 0.1, while the correspond- 
ing effective mass is shown in Fig. lb. Here we have suppressed any data 
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point which has error (jackknife with single elimination) larger than its 
mean value. Similarly, the eigenvalue A~(t) corresponding to the lowest 
J p = l/2~ state and its effective mass are plotted in Fig. 2a and Fig. 2b. 

First, we observe that the effective mass of the lowest J p = 1/2+ state 
in Fig. lb attains a plateau for t e [9, 14], and its mass can be extracted by 
single exponential fit (Fig. la). Similarly, the effective mass of the lowest 
J p = 1/2" state in Fig. 2b attains a plateau for t € [11, 17], and its mass 
can be extracted by single exponential fit (Fig. 2a). 




(a) (b) 

Figure 2. (a) The eigenvalue A~ (t) of the lowest negative parity state, for m u a = 0.1. 
The solid line is the single exponential fit for 11 < t < 17. (b) The effective mass of 
A~(t) in Fig. 2a. 

In Fig. 3a, the masses of the lowest lying J = (l/2) ± states are plotted 
versus mj, where all mass fits have x 2 /d-0-f- < !• Using the four smallest 
masses (i.e., with m u a = 0.03,0.035,0.04,0.045) for chiral extrapolation 
(linear in m \ ) to physical pion mass = 135 MeV, we obtain the mass 
of the lowest lying states: m(l/2~) = 1424(57) MeV, and m(l/2+) = 
1562(121) MeV, in agreement with the masses of rriK + mjv — 1430 MeV, 
and 9(1540) respectively. 

5. Distinguishing the KN scattering states 

Now the question is whether they are scattering states or bound states. In 
order to obtain the mass spectrum of KN scattering states, we consider the 
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time correlation function of KN operator without any exchange of quarks 
between K and N in its propagator, i.e., the interaction between K and N 
is only through the exchange of gluons. Explicitly, 

C KN (t) = (J2(N(x,t)N(d,0)) f (K(x,t)K(d,0)) f \ 

\ x I u 

where N = [u T C75d]d, and K = s-f 5 u. 

The masses of lowest lying KN scattering states are plotted in Fig. 
3b. For the J p = l/2~ state, using the four smallest masses for chiral 
extrapolation to m x = 135 MeV, we obtain m,KN0-/2~) = 1433(72) MeV, 
in agreement with the mass of vhk + rn^ — 1430 MeV, the KN s-wave. 
Further, its mass spectrum is almost identical to that of the lowest J p = 
l/2~ state of O in Fig. 3a, for the entire range of m u . Thus we identify the 
lowest J p = l/2~ state of <d(ududs) with the KN s-wave scattering state. 

On the other hand, for the J p = 1/2+ state, its mass is much higher 
than the naive estimate \Jm 2 K + (2tt/L) 2 + \Jm 2 N + (2-k/L) 2 , where L is 
the lattice size in spatial directions. This suggests that the KN p-wave (in 
the quenched approximation) in a finite torus is much more complicated 
than two free particles with momenta px = —Pn = 2nei/L. Further, the 
mass of KN p-wave scattering state in Fig. 3b is always larger than the 
mass of the J p = 1/2+ state in Fig. 3a. In particular, for m u < m s , the 
former is significantly larger than the latter. This seems to suggest that 
the lowest J p — l/2 + state of Q(ududs) is different from the KN p-wave 
scattering state. In other words, it is likely to be a bound state with mass 
1562(121) MeV. If it is identified with 9+(1540), then it predicts that the 
parity o/6 + (1540) is positive. 

6. Concluding remarks 

It is vital to re-confirm the above picture with a larger lattice, smaller lattice 
spacing, and higher statistics, especially for the regime m u < m s , which 
is crucial for the chiral extrapolation. Further, to ensure that the lowest 
lying J p = 1/2+ is not a scattering state, one can also compute its spectral 
weight for two different volumes 4 , since for a scattering state, its spectral 
weight is inversely proportional to the volume. Even if it is confirmed to 
be a bound state, one still has to find out whether its decay width could 
be as small as 20 MeV, compatible to that of 9 + (1540). Our present data 
only shows an unambiguous signal of pentaquark (ududs) resonance around 
1540 MeV, with S = +1 and I(J P ) = 0(1/2+). 
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Figure 3. (a) The masses of the lowest lying states of @(ududs). The solid lines are chi- 
ral extrapolation (linear in m^) using four smallest masses, (b) The masses of the lowest 
lying KN scattering states. The solid line (for J p = l/2~) is the chiral extrapolation 
using four smallest masses. 
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